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In most of the Nambu-Jona-Lasinio(NJL)-type models, realizing the hidden chiral symme- 
try, the existence of a scalar particle a is needed with a mass mo- = 2mq, as a partner of the 
Nambu-Goldstone boson tt. However, the results of many analyses on tttt phase-shift thus 
far made have been negative for its existence. 

In this paper we re-analyze the phase-shift, applying a new method, the interfering ampli- 
tude method, which treats the T-matrix directly and describes multi-resonances in confor- 
mity with the unitarity. As a result, the existence of a has been strongly suggested from the 
behavior of the tttt — > tttt phase shift between the tttt- and the KK- thresholds, with mass 
= 553.3 ± 0.5stMeV and width= 242.6 ± 1.23tMeV. The most crucial point in our analysis 
is the introduction of a negative background phase, possibly reflecting a "repulsive core" in 

TTTT interactions. 

The properties of /o(980) are also investigated from data including those over the KK 
threshold. Its mass is obtained as 993.2 ± 6.5st ± 6.9sj,sMeV. Its width is about a hundred 
MeV, although this depends largely on the treatment of the elasticity and the tttt KK 
phase shift, both of which may have large experimental uncertainties. 



§1. Introduction 

Chiral symmetry has a long history in particle physics and has played a central 
role in understanding the spectroscopy of hadrons with light flavors: In particular 
it is now widely accepted that the 7r-meson (or pseudo-scalar meson nonet) has the 
properties of a Nambu-Goldstone boson, corresponding to the dynamical breaking 
of chiral symmetry (D^SB) in the massless lirnit of QCD. 

In the Nambu-Jona-Lasinio (NJL) modelEP and its extended version (ENJL)o 
adapted to the quark model, which simply realizes the physical situation of D^SB, 
the existence of the cr-meson (or scalar meson nonet) is predicted as a chiral partner 
of TT (or a pseudo-scalar meson nonet), with iso-spin 1=0 and mass=2mg (rn-g being 
the constituent quark mass). 

It is also well-known that, in the investigation of nuclear forces with one-boson- 
exchange potential, eP an iso-singlet scalar particle with mass of about SOOMeV, which 



In the ENJL model is predicted the existence of NG boson, scalar meson, and moreover, 
vector and axial- vector meson nonets. 
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can possibly be identified with the a (or a strongly correlated two-pion states'), is 
necessary. However, the existence of o" as a resonant particle has not yet been 
generally accepted. A major reason for this is due to the analyses of vrvr uhase-shift 
obtained from the high-statistics data of a CERN-Miinich experimentEP in 1974, 
in which the phase shift 5q up to m,r7r = ISOOMeV turned out to be only 270°. 
Accordingly, after subtracting a rapid contribution of the resonance /o (980) (180°), 
there remains only 90°, which is insufficient for a around m=2mq=500 ~ 600MeV, 
and many analyses thus f ar m ade on the phase shift have yielded conclusions against 
the existence of (T.B''IZI''Q' hi 

Reflecting this situation, the non-linear realization of chiral symmetry now seems 
to be widely believed rather than the simple linear representation in the NJLmodel. 
In this case, the existence of a (or a scalar nonet in ENJL) is unnecessary, til which 
is represented as a function of tt (or a pseudo-scalar meson nonet). 

On th e o ther hand, the possible existence of a has been siigrested from various 
viewpoints US' both theoretically and phenomenologically.ll^^'tj^'llZP'lls* In particular 
the importance of a in relation with the DxSB has been argued extensively by 
Refs|l|),|g),^. 

Here, it may be worthwhile to note that in a recent pp central collision experi- 
ment, an event concentration in the 1=0, 5"- wave vrvr channel is seenE3) in the region 
of m^rTr around 500 ~ 600MeV, which is too large to be explained as a simple "back- 
ground" and seems strongly to suggest the existence of a (a brief review is given 
in §4(B)). 

Considering all the above situations, a more rigorous re-analysis of yrvr-phase 
shift seems to be necessary. 

In this paper, the authors will perform this by applying a new method of analysis, 
the Interfering- Amplitude (I A) method, which is able to treat partial- wave multi- 
resonances directly in the T-matrix (instead of in the most conventionally used /C- 
matrix) in conformity with unitarity in the final state vrvr-inter actions. As a result, 
evidence for the existence of the a-meson is shown from the characteristic behavior 
of the vrvr — > vrvr phase shift under the KK threshold. A keystone for this is the 
introduction of a negative background phase, due to an unknown repulsive force 
between pions, which might suggest the existence of a repulsive core in the vrvr 
interaction. Here it is notable that the existence of a similar negative background 
phase has been established phenomenologically in the case of the A^-A^-interaction. 

In §2, the lA Method is introduced. Its relation to other approaches, such as 
the /C-matrix method, is discussed in §4(A). 

In §3, the tttt scattering phase shift is re-analyzed by applying the lA method. 
It is especially shown that the characteristic features of phase shifts from the vrvr- to 
the KK- thresholds, including a small bump around 750MeV, are reproduced well 

**■' However, see the analyses which suggest the existence of a. 

*' See, e.g., Refs. ^ and|l3|). 

**•' Owing to Watson's final state interaction theorem, the nn production ampfitude of all types 
of reactions is closely related to the amplitude of tttt scattering. Accordingly, if this concentration 
is actually due to the cr-particle production, the corresponding phase shift ought to be observed in 
TTTT scattering. 
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with our new background phase and the a with mass of 555MeV. 

The properties of /o(980) are also examined with the data including those jjver 
the KK threshold in this section. Its width attains a value from tens MeV™ to 
hundreds of MeV,lly^ largely depending upon the treatment of the elasticity and the 
vrvr KK phase shift, both of which have large experimental uncertainties. 

In §4(C) we give some arguments on the repulsive core in the vrTr-interactions, 
a possible physical origin of the negative-phase background. It is pointed out that 
our best fit value of the core radius is almost equal to the structural size of the pion, 
showing a similarity to the case of the A^-A^ interaction, where the core radius has 
an order of the size of a nucleon. 

In §4(D) the existence of new types of 0^+ and 1^+ nonets is suggested. These 
are outside of the normal L-excited meson nonets, a seems to have properties such 
that it should be assigned as a member of this new type of scalar nonet. 

§2. Partial-wave multi-resonance and interfering-amplitude method 

In extracting resonances due to physical particles from scattering amplitudes, we 
usually fit partial wave amplitudes with the Breit-Wigner(BW) form. The scattering 
matrix S must satisfy unitarity, which means conservation of probability from initial 
to final states: 

SS^ = S^S = 1. (2-1) 

We shall treat the vrvr 5- wave scattering in case of two relevant channels: vrvr and 
KK, denoted as channels 1 and 2, respectively. In this case the p matrix has a 
diagonal form 

p = diag{- ^^^^ ^^^^ ^ 



'8vry^' Svr-^/s ' 



IpiI = \/ J - ml, |p2| = W J - m\, (2-2) 



The partial wave scattering matrix is defined as 



5ji = {2^f5^'\P, - P^)J^f: ^P/(n • n04, (2-3) 

V P:iPi 1=0 ^ 

where the suffix j(i) represents the final (initial) channel, n'(n) is a unit vector of 
3-momentum of the first particle in the final (initial) system. The bold suffices j, i 
include momenta of the system. 

The unitarity equation (2T) is represented, in terms of 5"', as 

S\Jl=hu / = 0,1,2,.... (2-4) 
The 5' matrix is symmetric due to time reversal invariance 
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The T-matrix and its partial wave component are usually defined, for relevant 
systems of two spinless particles as 

5ji = 5ji-i(2^)V^)(P,-P,)'?ji, (2-6a) 
Tji = ^(2/ + l)P,(n.n0'7;',. (2-6b) 



The is related to T' as 



S'j, = dji - 2i^J-pir^i. (2-7) 



Equations (2-4) and (2-5) are equivalent, respectively, to the following conditions on 

i{Tl - T^) = 2(rVr't),,, (2-8) 

Vi = '^r (2-9) 

In the case with only one resonance, these conditions are satisfied by the BW ampli- 
tude. However, in the case of "multi-resonance scattering" (when many resonances 
exist with the same quantum number), if we represent the partial wave T matrix as 
a simple sum of corresponding BW formulae, they are not satisfied. 

Here wc give a new method for describing the scattering amplitude in conformity 
with the unitarity relation in this case. 

We shall treat directly the partial wave amplitude defined as 

4i = -VPjP^'^ji^ (2-10) 

which is related to 5' as 

S'j, = 6ji + 2ia'j,. (2-11) 

In the following we treat only a^j (or S'jj) with I = 0, and the superscript I is 
omitted. 

2.1. one-channel multi-resonance scattering amplitude 
First let us consider the one channel case. 

A basic idea of our method is that a phase shift 6{s) of the scattering amplitude 
is essentially only due to physical particle resonances through the BW formulae, and 
the phases of multi-resonances contribute additively to S{s). Thus it is represented 
as 

o (1) (2) (3) 

S{s) = S'^ = S + S + S +■■■ , (2-12) 

where 5 is a real phase coming from the i-th resonance. 

In the case with one-resonance, the scattering matrix is given as 

(1) (1) (i) 

S=e^'^ = l + 2i a, (2-13) 



An Analysis of inr- Scattering Phase Shift 
by the BW formula 

(1) -VsrUs) , , IpiI «^ , , w 



2 



a 



, VsrUs) = 9 1= piis) g 1 . (2-14) 



(1) 

The 6 is real, as is seen by substituting Eq.(2-14) into Eq.(2-13), and the uni- 
tarity relation SS* = 1, Eq.(2-4), is easily shown to be satisfied. 
In the two-resonance case the scattering matrix, 

S = e'^'^'-'^ = l + 2ia{s), (2-15) 

is given as a multiplicative form of the "respective resonance S-matrices" S= e^^^ = 
1 + 2i a as 

(1)(2) 

S=SS. (2-16) 
The unitarity relation SS* = 1 is easily seen to be satisfied by the "unitarity of 

respective S*-" The scattering amplitude a{s) (= —pi{s)T{s)) is represented, in 
terms of the respective BW resonances, as 

a{s) =a {s)+ (s) + 2i (s) (s). (2-17) 

This amplitude consists of a simple sum of respective BW resonances and their cross 

term [2i ^a! (s) (s)). The latter, which looks like an "interference" of the former 
two resonances, appears in the amplitude. Thus we call this method of constructing 
the amplitudes the Interfering-Amplitude (lA) method. 

Extension to the three-resonance case is straightforward, and the formulae are 
given as 

(1) (2) (3) (l)(2)(3) 

5=e2^(^+^+^)=yyy, (2-i8) 

(1) (2) (3) (1)(2) (2)(3) (3)(1) (1)(2)(3) 

a[s) = a + a + a +2i[ aa + aa + aa) — Aaaa. (2-19j 

A background phase is added to the phase shift Eq.(2-12), if necessary: 

S = 6^^ + S^^. (2-20) 

The unitarity relation is satisfied for any real 6^'~^ . In this case the scattering am- 
plitude is given by a background amplitude a^*^ and the resonance amplitude a^^ 
(Eqs.(2-17) and (2-19)) as 



a = a 



BG + ^Re ^ 2m^^a«^ = a^^ + a^^e^^^'''' , (2-21) 



*^ Here we apply a relativistic form with a modification, taking into consideration the indefi- 
niteness of the masses of unstable particles. Decay widths of resonances will be estimated using this 
form in 53. 
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where 



2i 

Application to the greater than three-resonance case can be done in a similar way. 

2.2. two-channel multi-resonance scattering amplitude 

Next we construct an explicit form of the scattering amplitude in the two-channel 
case. We treat scattering matrix elements ^n, 5i2 and S'22. The element 521 is equal 
to by time reversal invariance (Eq.(2-5)). On these three elements is imposed 
the unitarity relation (Eq.(2-4)) 

I'S'llP + |S'i2p = |522p + |'S'l2|^ = 1, 

SnSi2 + S12S22 = 0- (2-22) 
The form of Sij, satisfying the above equations, is written as 

511 = r? e^^^ 

512 = ^l-r?V('^+i), 

822 = Tje'"'^', (2-23) 

in terms of the phase shift S{S') in the tv'!t{KK) scattering, the corresponding elas- 
ticity r]{= 77'), and the phase (f) in the tttt — KK scattering. These phases are 
constrained by 

6 + 6' = (2-24) 

In the one-resonance case, the scattering matrix elements satisfying the unitarity 
relation Eq.(2-22) are given in terms of the BW amplitude as 

(1) (1) (1) (1) (1) (1) 

5ii = l + 2zaii, 5i2=2zai2, 522=l + 2i a22, (2-25) 



Qj in 



(1) -v'v^^('i)v^'"a) 



Phases 6 and (p are obtained by comparing Eqs.(2-25) and (2-26) with their definition 
Eq.(2-23), and the elasticity 77 is given as 



(1) 
V = 



1 , ^/^^(\)(^)V^^fi)(^) ,,27) 



In the two-resonance case, extending the method in Eq.(2-16) to this case, the 
diagonal elements of the scattering matrix, 5ii and ^22, are represented in a multi- 
plicative form of respective resonance 5-matrices 

(1) (2) (1) (2) 

Su — 5ii5ii, S22—S22S22, (2-28a) 
'L = %'e-'Vl2.4'e-', (2.28b) 
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where the elasticity i] in each diagonal channel is represented in the multiplicative 

(1)(2) (1) (2) 

form ij = 11 V , and the phase shift as a sum 6 = 5 + S , following from its physical 
meaning. Corresponding amplitudes are given in a form similar to Eq.(2-17) in the 
one-channel case: 

Su-l (1) (2) (1) (2) 

an = — — — = a 11 + an +2t ana n, 

S22-I (1) ^(2) ^^.(1) (2) 

022 = — ^. — = o 22 + o 22 +2i a 22 o 22 ■ (2-29) 

The non-diagonal S'-matrix element is determined through relations (2-23) and (2-24) 
with the diagonal elements given above. First, the phase shift (f) of S12 is equal to 

(1) (2) 

a sum of respective contributions of resonances: (j) =<p + • So it becomes of the 

form 5i2 = 2i ^a 12 x ^au x(real function). The absolute value of 5i2 is determined 
by Eq.(2-22) 

„ .(1) (2) VI -J?^ o- fnon\ 
S12 = -i a 12 o 12 -(I)— = 2zai2. (2-30) 

I o 12 o 12 I 

The unitarity Eq.(2-22) is guaranteed by the above procedure to derive S'12. 

Extension to the three-resonance case can be done in the same manner. The 
results arc given as 

_ (1) (2) (3) 

S ji — S a S a S a tor t — 1)2, 



.(1) (2) (3) y^l^ 



St2 = — 1 O 10 O 12 O 10 , , , , , , , 

1^ iZ IZ LA (3) ' 

I O 12 O 12 O 12 I 

{1)(2)(3) 

r]= 11 11 11 , (2-31) 



_(1) (2) (3) {I) (2) (2) (3) (3) (1) . (1) (2) (3) 

Oil — O ii + O jj + Oj ii 'T^^iy CI jj O jj -|- O jj O jj -f- a a j^^j 4 O j j O i j O i j , 

012 = (2-32) 

Here, it may be worthwhile to note that all phases, 6[s),5'{s) and 0(s), are given as 
a sum of contributions of respective resonances as in Eq.(2T2). Application to the 
greater than three-resonance case can be done in a similar way. 

Background phases are added to phase shifts if they are necessary. In the two- 
channel case, in order to satisfy time reversal invariance, that is S'12 = S'21, back- 
ground phases are introduced as 



-<BG 




(2-33) 
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where S^'^ {5'^^) denotes a phase in the vrvr {KK) channel. This leads to the 
following change of the S-matrix phase: 

oRe QRe2iSBG 

^11 ~^ '-'11 6 

cRe _^ cReJiSBG+s'BG) 
'-'12 '-'12 6 

S^f -5te2^^'"". (2-34) 

In the actual analysis we introduce a background phase only in the vrvr channel, that 
is, 6'^^ = 0. In this case we obtain 

ail — flu + aiie , 
ai2 — ai2 e , 

,Re 



022 = a^l, (2-35) 



where af/^ — • 

Finally, in this subsection we rewrite our formulae of the two-channel scattering 
matrix in a simpler form in terms of 

n=s- Mg) + + V~srl^{s)). (2-36) 

In the one-resonance case they are given as 



(1) (1) (1)* (1) / (1) (1) (1) 

5ii=W/ d, S22=n I d, Su = -i\l\ d P-l / d • (2-37) 

In the multi-resonance case they are 

s„=n»/'2.&=nS'/2. 



'12 



The above form of the one resonance scattering matrix coincides with the one given 
in Ref. ^), whose method can be naturally extended to the multi-resonance case and 
gives the same form as Eq.(2-38). 

§3. Analysis of vrvr-scattering phase shift 

Figure 1 shows the 5-wave, 1=0 vrvr phase shift reported in the analysis of the 
CERN-Miinich data, & which is widely accepted as a standard of the phase shift up 
to 1900MeV. 

It shows a rapid step up by 180° near the KK threshold, following a slow increas- 
ing which reaches 90° slightly below DOOMeV. Above the KK threshold, it continues 
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to grow slowly up to about 270° at 
ISOOMeV and 360° at 1600MeV. The in- 
crease in the phase shift by 180° corre- 
sponds to one full-resonance. Ac- 
cordingly, the phase shift up to ISOOMeV 
(~ 270°) has been interpreted so far as 
owing to the existence of two resonances, 
i.e. a narrow /o (980) (full-contribution by 
180° because of its narrowness) and a very 
broad /o(1300) (half-contribution by 90°). 
The latter, which corresponds to e(1200) 
in 1976 PDG©, replaced the pole pre- 
viously listed below 600MeV, and a has 
been omitted from the PDG table since 
the early 1970's. 

However, we have found some diffi- 
culties in the conventional interpretation 
mentioned above. We have made a fitting 
of phase shift below ISOOMeV following 
this interpretation, for which the results 
are given in Fig.l. As is clearly seen, the 
significant structure of a small bump of 
the phase shift around 750MeV cannot be reproduced by this interpretation. This 
structure clearly requires certain sources around this mass region. Moreover, it seems 
that the behavior of phase shift below 600MeV, which will be given in §3.1, is not 
reproduced by this fitting. 

To solve these two problems, in analogy with the case of nuclear forces, we have 
introduced a negative background source, which works to pull back the phase shift. 




500 



1500 



Fig. 



1000 

M&P&P (MeV) 

1. TTTT scattering phase shift. The fitting 
by the conventional interpretation with 
/o(980) and /o(1300) is also shown. 



^c|-f*7r7r I ; 



(3-1) 



where rc is a constant parameter with dimension of length, and Fj^n is equivalent 
to pi in Eq.(2-2). Then the phase and amplitude are described, including the back- 
ground, as in Eq.(2-20) and Eqs.(2-21) and (2-S5), respectively. The source of this 
negative phase background must be a "repulsive" force between pions. We will give 
some arguments on its possible physical origin in §4(C). 

In the following, we perform fittings on the various data using the lA-method 
with the background phase. In §S.l, an analysis of the vrvr phase shift under the 
KK threshold is presented, focusing on the possible existence of (t(555). In §3.2, 
the data including those over the KK threshold, the vrvr phase shift, the elasticity 
and the vrvr — > KK phase shift, are examined all together in order to investigate the 
properties of /o(980). Finally, in §3.3, a more detailed study is made on the shape 
of the negative background, trying a fitting in the entire relevant mass region. 



*-* We assume here that our relevant resonance particles couple with the tvtv channel strongly 
enough to make this phase shift. 
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threshold — Existence of the cj(555) particle and its 
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3.1. Fitting 6 below the KK 
properties 

We first perform a fitting for the 
TTTT — s- TTTT phase shift ((5) below the KK 
threshold (280MeV-980MeV). In this re- 
gion, we suppose the existence of two res- 
onances, a and /o(980), and introduce the 
negative background mentioned above. 

We apply the one-channel two- 
resonance formulaO given in Eqs.(2-14), 
(2-15) and (2-17) with the background 
phase described in Eqs.(2-21) and (3-1). 
Necessary parameters to be determined 
are the following five: Mc(= 

(1) (2) 

9c7r7r(=5i), Mg{= M(2)), gg-K-K{=9i) and 
fc] 9-K-K is the coupling constant to tttt 
channel, where suffix c and g represent a 
and /o(980), respectively. Tc is the "ra- 
dius" of the negative background. Cou- 
pling to KK is not needed, because 
we only treat the phase shift below the 
threshold here. 

The CERN-Miinich analysis provides 
data in the mass region over GOOMeV. 
For acquiring the property of o", we sup- 
plement the data from the vrvr threshold 
(2S0MeV) to eOOMeV by Srinivasan et al. 
c§) and by Rosselet et alx3 

Figure 2 shows the result of our fit- 
ting. Values of resonance parameters ob- 
tained are listed in Table I. The jmpst notable result is identification of cr with the 
low mass of 553.3 it 0.5,/MeV. 



-90 



1 1 1 1 1 1 1 
+ G.Grayeretal. 




♦ L.Rosseletetai. 
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Mspsp (MeV) 



1200 



Fig. 



2. Fitting of 5 below the KK thresh- 
old (shown by an arrow) with two reso- 
nances. Behavior of the negative back- 
ground phase is also shown (dashed line). 
The parameters used are described in the 
main text. Their best fitted values are 
listed in Table I. 



According to the result of another fitting to 
data including those over the KK threshold, which will be presented in the next 
sub-section, the KK-coupling of C7(555) turned out to be negligible. Thus the total 
width 242.6 lb 1.2siMeV owes only to its vrvr coupling. The fine structure of the phase 
shift around 280 ~ DOOMeV, especially the bump around 750MeV, is reproduced 
quite well. It is due to the cooperation of the phase increase by (t(555) and decrease 
by the negative background. The radius of the negative background turned out to 
be 3.46GeV-^ 



*' See the remarks to be given in §3.4. 
**' There are several choices of the relevant relativistic Breit Wigner forms. The mass value is 
dependent to some extent upon these choices. 
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Table I. Parameters in the best fit below the KK threshold with 2 resonances and the negative 
background phase. Properties of /o(980) in this table may not be definitive, due to omitting data 
over the KK threshold. For the property of /o(980), see §3.2 and Table II. The value of ff^-jj- is 
quoted from the result of fitting data including those over the KK threshold. The decay width 
r and the "peak width" r'"' are defined, respectively, as n = J^'^ dsri{s)s^^'^ r^'°\s) /{tvKs - 

Mff + sr!°\sf]) and I^'"' = Fii^s = Mf). The relation of g and r(s) is given in our BW 
formula Eq.(2. 14). If we use another form of BW formula, instead of (2.14), as -Mr{M^)/{s~ 
+iMr{M'^)), the pole position is simply given as Spoie/M—M — ir{M'^). See the discussions 
in §4(A). 





mass(MeV) 




9kkI 3'^"" 




r^f'^MeV) 


a 


553.3±0.5 


3336±12 


0. 


242.6±1.2 


349.3±2.5 


(/o(980)) 


(970.7±2.2) 


(1768±24) 







2 BW(GeV-^) 




3.46±0.01 



3.2. Fitting of 5, tj and (j) over the KK threshold — properties of fQ{980) 

In the following we will investigate the vrvr vrvr phase shift (6), the elastic- 
ity [r]), and the vrvr KK phase shift {(j)) over the KK threshold all together, to 
study the properties of f n (980). We shall try to fit data in the mass region from 
eOOMeV to ISOOMeY H with three resonancesLj2H-,/o(980) and /o(1300). We apply 



the two-channel three-resonance formula given as Eqs. (2-31) and (2-32) 

with a background phase (2-35) and (3-1). Thus the total number of parameters 

(1) (2) (3) 

is 10: Me(= M(i)), Mg{= M(2)), Mo(= M(3)), gc^^(=5i), <?9^^(=ffi), go^^(=9i), 

(1) (2) (3) 

9cKKi=3 2), 9gKKi=3 2), 9oKKi=3^) ^nd r^. 

Data treatment 

Our treatment of data is as follows: 

(a) (5 over eOOMeV bJ is quoted from the CERN-Miinich analysis, i As can be seen 
in Fig.l, one data point just on the cliff of /o(980) exhibits strange behavior. It 
looks like a sharp peak, which might not happen in a physical phase shift. Including 
or excluding this point seriously influences the results on /o(980). Thus we try the 
following two types of fitting, i.e., one which does not include this data point (fit A) 
and another which uses all data points (fit B). 

(b) ?7 is quoted from the same analysis. Also there is a point corresponding to the 
point referred in (a) just in the narrow dip about IGeV, which seems to require quite 
a large inelasticity. In fit A we will not use this point and in fit B we will use it. 



***•* See the caption of Fig. 3. 

****' Below llOOMeV, the intermediate states in the final state interaction may be limited to the 
TTTv and KK channels. We will treat only these two channels up to 1300MeV. 
,****) discussion to be given in §3.4. 

*' The procedure omitting the data below 600MeV has some influence on the mass and width 
of a. It implies that the behavior of the background phase is not as simple as in Eq.(3T). This 
problem will be discussed in §3.3. Anyway, as far as inquiring into the physical problem over and 
around the KK threshold, the low energy region below GOOMeV is not so important. Thus we use 
only the CERN-Miinich data for fitting here. 
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Fig. 3. The fitting for (a) S, (b) rj, and (c) over tfie KK tliresfiold witfi tliree resonances (fit A: 
witliout tlie two data points (mentioned in tlie text), fit B: with all data points, fit C: omitting (p 
behavior). Data in 600MeV ~ ISOOMeV shown by arrow are used for the fitting (the data over 
ISOOMeV are omitted because we must consider the possible effect due to the other resonances 
in that mass region). 



(c) (p is quoted from Ref.|^. As can be seen in Fig. 3(c), different conflicting results 
are reported on the ^-behavior. Thus we shall n£rform a compromised fitting: we 
use a simple mean of the values by Etkin et a/.E2P(<S) and by Cohen et a/.EiP(C), a 
square-mean of their errors as statistical errors, and their differences as systematic 
errors : mean = {£ + C)/2, error = y/ (AS)''^ + {^C)'-^ + \£ — C\. 
As will be shown later, our result reproduces that of Etkin et al. rather than that 
of Cohen et al. The fitting by Morgan and PenningtonQ' also exhibits a similar ten- 
dency. Our result does not change significantly when using only the data of Etkin 
et al. 

Phase 6 supplied from negative background 

Fitting of 6 is shown in Fig. 3(a). Concerning our relevant problem, there is no 
significant difference between the cases in which we use fit A and B. 

Note that we have enough space of S for 3 resonances: a, /o(980) and /o(1300) 
below 1300MeV, due to the contribution of the negative background. The other 
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Table II. Parameters in the best fit over the KK threshold with 3 resonances and the negative 
background phase. The results on a in this table are somewhat different from those in §3.1, due 
to omitting the data below 600MeV. This problem is solved in §3.3. For reliable results on a, 
see Table I. The results on /o(1300) may not be definitive due to the uncertainty in the choice 
of background phase. See the main text. On the definition of F and r^^\ see the caption in 
Table I. 





mass(MeV) 


p^^(MeV) 


9kkI 9'^'^ 




(598.6±15.6-29.6) 


(4126±348) 


(0.0±0.349) 


/o(980) 


fit A 


993.2±6.5±6.9 


1680±91 


1.48±0.23 




fit B 


982.3±2.2±2.4 


1442±82 


2.82±0.22 




fit C 


987 


1947 


3.77 


(/o(1300)) 


(1313±36) 


(5185±272) 


(0.286±0.079) 





rtot(MeV) 


r r — 




(325.1±36.3) 


(325.1±36.0) (0.0±2.8) 


/o(980) fit A 
fit B 
fit C 


78.3±11.5 
114.5±15.5 
344 


55.70±6.1 22.7±9.6 
46.1±4.9 68.4±14.6 
73 272 


(/o(1300)) 


(372.0±34.5) 


(355.9±34.2) (16.1±9.9) 




/^r(MeV) 




(^) 


(505.0±89.3) 


(505.0±89.3) (-) 


/o(980) fit A 
fit B 
fit C 


67.9±9.4 
40.5±4.7 
73.5 


54.4±6.1 13.5±7.4 
40.5±4.7 
73.5 


(/o(1300)) 


(420.6±46.2) 


(398.6±44.1) (22.0±13.8) 





3 BW(GeV-') 




3.76±0.27 



resonances, cl such as a glueball candidate can obviously exist with the mass around 
1500 ~ leOOMeV. 

Results on /o(980) 

The mass of /o(980) is obtained as 993.2 ± 6.5st ± 6.9sj^sMeV(fit A), where the 
systematic error is estimated as a probable uncertainty {2a) of Tc (see §3.3). 

On the other hand, its coupling to KK channel varies a great deal, depending 
on whether we adopt fit A or B (Table II). 

In Fig. 3(a) the 6 fitting for fit B is also shown. As can be seen, the fitting over 
IGeV is worse than that for fit A, due to the infiuence of the data point on the cliff, 
referred to above. 

The elasticity plots for fit A and B are shown in Fig. 3(b). If we take the point on 
the cliff into account (fit B), the fitted curve of the elasticity passes through the lower 
region, and the (/)-fitting becomes slightly worse (Fig. 3(c)). Thus the KK width and 
the total width in fit B must be larger than those in case of fit A, respectively. If 
we intend to reproduce that point exactly, the KK width and total width become 
272MeV and 344MeV, respectively (fit C). They seem to be consistent with recent 

*' However, they may be "partial resonances", which couple with the tttt states only partially, 
when the corresponding phase goes up and then down, and finally giving no contribution to the 
phase shift. 
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values obtained by Zou et al.^ But in this case it is difficult to reproduce the cf) 
behavior as shown in the figure. The obtained width of /o(980) without this point 
(fit A) is 78.3 ± ll.SgtMeV, which is consistent with a value given by Morgan and 
Pennington. lP 

It should be noted that the results on properties of /o(1300) in Table II are not 
definitive, because they strongly depend on how to choose the background phase, 
and also on the properties of other possible /o with higher masses. 

3.3. Fitting for the entire relevant mass region — possible form of negative back- 
ground 

As was mentioned in the previous two sub-sections, some differences of the ob- 
tained results between the fittings in the different mass regions suggest that the 
shape of background phase is not as simple as Eq.(3T). 

From a consideration on the origin of 
negative background, considering an anal- 
ogy with the case of nucleon-nucleon inter- 
actions (to be given in §4(C)), we consider 
it reasonable that the radius of the nega- 
tive background Tc becomes smaller with 
the increase of the pion momentum. Thus 
we shall further try another type of fitting 
through the whole relevant mass region to 
determine a possible shape of the back- 
ground: _ 
(a) We will introduce an exponentially de- £ 
creasing factor to the radius of the back- S 
ground Tc (Eq.(3T)), 



360 



270 



180 



rcexp{-b\Pj, 



(3-2) 



where the parameter b also has dimensions 
of length. 

(b) Parameters giving the properties of 
three resonances are fixed, i.e., those for 
a are taken from the values obtained from 
fitting below the KK threshold, while 
those for /o(980) and /o(1300) are from 
the values over the threshold. 

The result of fitting is shown in Fig. 4, 
where the best fit parameters of the back- 
ground core are given in Table III. The 
fitting looks quite good through the en- 
tire relevant mass region. 
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4. Fitting for S through entire mass re- 
gion with "soft" repulsive core. The pa- 
rameters used for a are quoted from Table 
I, while those for /o(980) and /o(1300) are 
from Table II (fit A). 
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Table III. Parameters in the best fit for tlie core of Eq.(3.2). 



rc(GeV-') 


b(GeV"') 


5.36±0.01 


0.480±0.009 



3.4. Remarks on the method of our analysis 

In performing the analysis in this section we have actually applied an analytically 
discontinued (A.D.) density matrix, whose elements have non-zero values only over 
the threshold of the corresponding channel as 

/ lPl| IP2I N 

p = diag{- — — -), 



\pi\ = ^l-ml9is-4mi), = ^ - - m^eis - 4ml) (3-3) 

(correspondingly in fitting of §3.1 below the KK threshold we have applied the one- 
channel formula, and in fitting of §3.2 applied the two-channel formula, since only 
the |pi| remains non-zero for the region below the KK threshold). 

Thus the resonance parameters obtained in this analysis are considered to repre- 
sent the particle properties in the zero-th approximation, which are determined only 
through the quark-gluon color-confining dynamics, omitting the "residual" strong 
interactions among color-singlet hadrons. However, this treatment violates the "sa- 
cred" requirement of analyticity of the scattering amplitude. Applying an analyti- 
cally continued (A.C.) density matrix is considered to take into account the effects 
of virtual tttt and KK channels in extracting the resonant-particle properties. 

We have estimated these effects, actually applying the A.C. density matrix: The 
result of re-analysis of data (including all of 6, rj, and (j)) in the region 600<y^<1300MeV 
produces little change in the properties of (t(555) and /o(980). This seems to be due 
to the zero coupling of a to the KK channel {of. Table I) and the small width of 
/o- Then, in order to estimate the effects on the cr-properties by applying the A.C. 
density matrix for /o, we have performed fitting under the KK threshold by tak- 
ing the /o properties as obtained above and by setting gp,-^{a)=0. The results are 
M<^=568.9±0.5MeV and r^=250.5±3.7MeV. The smah differences between these 
values and those in the case with A.D. density matrix (Table I) owe to the A.C. 
effect on /o, which has non-zero coupling to KK. 

In the analysis of this section we have not taken into account the Adler-zero 
conditioned) for the amplitude: The possible effects of this requirement on the prop- 
erties of a and of /o are considered to be small, although it may affect the values of 
parameters relevant to the energy region near the vrvr-threshold. 

Here we should point out that it is generally very difficult to take quantita- 
tively the effects of strong interaction into account, since hadrons are interacting 
"democratically" among various possible channels (remember the long history of 
strong interactions before QCD). Also we feel that the physical backgrounds of the 
analyticity of the scattering amplitude in the level of point-like hadrons should be re- 
examined presently, when the "global" properties of hadrons are determined through 
the confinement dynamics of the quark-gluon world and the hadrons have space-time 
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extensions. 

§4. Physical implications of our results and related problems 

(A) Relation of the lA method to the other approaches 

The S'-matrix must satisfy the unitarity (Eq.(2-4)) and be symmetric due to time 
reversal invariance (Eq.(2-5)). The corresponding conditions (2-8)^nd (2-9) are im- 
posed on the T matrix, and there have appeared various methods^ to parametrize 
the general solutions satisfying these conditions^ 

In the most conventional /C matrix methodslIP'EP the T matrix is first represented 
by the /C matrix as 

and is regarded as an analytic function on the Riemann sheets of complex s-variable, 
except for the right and left hand singularities and the poles. Then the T matrix 
for any real symmetric IC matrix automatically satisfies the unitarity and the time 
reversal invariance. 

Our relevant problem is to find pole positions of the T matrix, which are related 
to the masses and the decay widths of resonant particles. 

However, the parameters used in the IC matrix are not directly related to these 
properties of physical particles, and there seem to be many unpleasant features: 
Even the number of T matrix poles is determined by the choice of the /C matrix 
parametrization. Moreover, the contributions to the T matrix of the resonances and 
the non-resonant background cannot be separated in the simple /C matrix represen- 



tation. In the Dalitz-Tuan representation eIP of the T matrix, which is obtained by 
modifying the IC matrix representation, the background and resonance contributions 
are clearly separated, and the parametrization used in Ref.^) is equivalent to our 
Eq.(2-35) in the one-resonance case. 

On the other hand, in the w-method which was proposed by KatoE3) and by 
Fujii and Fukugita, & the pole positions of the S'-matrix are directly parametrized 
on complex Riemann sheets. They use the Le-Counter Newton equation for the S 
matrix in two-channel case, 

= dW(s)/d(s) for i = 1,2, 
SuS22-Sf^ = S''^\s)/d{s). (4-2) 

Here d(s), defined on the physical Riemann sheet /, includes the right hand singu- 
larities, while the d^"^^(s), d*^^-*(s) and d^^''^\s) are the analytic continuations of d(s) 
in the sheet II, IV and III, obtained by the change of signs of momenta pi — > —pi, 
P2 — *■ —P2 aiidi^3i,p2) — > {—Piy —Pi)- The four Riemann sheets are mapped onto the 
single w-planeE^ 

^= ; / ={Pl+P2). (4-3) 
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The poles of the S matrix correspond to the zeros of d(s) in the w-plane. The Breit 
Wigner resonance formula is obtained with d(s) 

2 2 

^ m^- m^ _ ^ _ ^ 

where the four zeros {uJr, and their inversion for the imaginary axis — w*, —a;*) 

satisfy certain constraints from the unitarity. Then, d(s) becomes equal to d (s) in 
Eq.(2-36), if we replace y/sri{s) = Vs-^gf by mr-{s) = m ^^^^ gf , and the final 
form of the S matrix in Eq.(4-2) coincides with our Eq.(2-37). 

In our case with ■y/sl^j(s) the 5-matrix has six poles. Since the ^-matrix in 
the lA-method is described by the simple Breit Wigner resonance formulae and is 
parametrized directly in terms of the masses and coupling constants, it seems to 
us not necessary to determine the pole positions of S'-matrix which have no direct 
physical meanings. 



(B) Central production of a (555) 



According to Watson's final state interaction theorem, the phase of production 
amplitude of all types of tttt production processes is determined through the tttt- 
scattering amplitude. Accordingly, if the existence of (t(555) is really observed in 
the vrvr-scattering, the corresponding signal must also be observed in any production 
processes. 

In the observations of the tttt production in pp central collision: pp — > pf{TT-K)ps 
with hundreds of GeV/c momentum incident proton bearri. a huge event concentra- 
tion around m-TTTr = 0.3 ~ IGeV has been clearly seen.iB''E3^'ii^ In the conventional 
treatment thus far made, this was regarded as a simple "background" . 

However, we findE3) that the /olQSO) peak slides down around 900MeV, probably 
due to "destructive" interference l1 with the concentration. We have performed a 
fitting with two Breit- Wigner resonances (/o(980) and another of lower mass and 
wide width), and found that "destructive" interference occurs with a relative phase 
between the two resonances, which seems to be consistent with the above theorem. 
This fact supports our interpretation of this concentration being due to production of 
cj(555). However, it is to be noted that the mass and the width obtained in Ref.22) 
may not be definitive, and the above result should be checked further. An exact 
analysis on the pp central collision will be reported soon.© 

Moreover, by comparing the vrvr-spectra in the pp central collisions and in the 
peripheral 7r~p charge-exchange processes, it is seen that the production probability 
of a"(555) in the former may be larger than that in the latter. This fact seems to 
be explicable by supposing that o"(555) is a member of the new type scalar nonet, 
which will be discussed in (D) below. 



*' In Ref.g8|) an interference of the /o(980) with the S'-wave background was introduced. In 
Ref.p^) it is first pointed out that this "interfering background" may originate from a physical 
particle. 
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Table IV. The value of hard core radius of tttt interactions, in comparison with the cases of a-a 
and N-N interactions. Structural sizes (EM r.m.s radii) of relevant particles are also shown. 
* value from Table II. ** value of soft core from Table III. A permissible minimum of the tttv 
core radius turned to be 0.5fm (2.5GeV~^) (see the text). 



system 




V< r'^ > 


a — a 
N-N 


2/m 
(lOGeV"') 
0.4 ~ 0.5fm 
(2 ~ 2.5GeV-i) 


1.6fm 
0.8fm 


TV — TV 


0.68fm 
(3.46GeV-^)* 

1.06fm 
(5.36GeV'^)" 


0.7fm 



(C) "Repulsive core" in tit: -interactions 

In our analysis of the vrvr-scattering phase shift of §3, leading to strong evidence 
for the existence of ct(555), introduction of a negative-phase background Eq.(3-1) has 
played a crucial role. We consider that a possible physical origin of this background 
is of a very fundamental nature and that a corresponding repulsive force may be 
reflecting a Ell Fermi-statistics nature of constituent quarks and anti-quarks of pions 
in the vrTT-system: A repulsive force due to a similar mechanism has been known 
since long ago to exist in nucleus-nucleus interactions. In the a-a system it isphe- 
nomenologically known that there exists a strong short-range repulsion0)'0'@'El* 
which gives a negative scattering phase shift of the same form as Eq.(3-1) with a cor- 
responding size of core r"" = 2fm (, while the r.m.s. radius, a "structural size", oLq 
is \J <r^ >a = 1.6fm). The physical origin of this repulsion was explained to he'^ 
due to the Fermi-statistics property of constituent nucl eon s of each a-nucleus in the 
a-a system. Then it may be reasonable to imagine, o^'cS'H by replacing the con- 
stituent nucleons in the a-particle with constituent quarks in the nucleon, that there 
may exist a strong short-range repu 
system. Actually it has been shown 



i n n uclear forces between nucleon-nucleon 



on ir 

P'cZP'E^P phenomenologically that there surely 
exists a strong repulsive core in ^Sq states. Extensive investigationsll^) on nuclear 
forces have been done by many physicists, and the size of the corresponding hard 
core radius was reported to be r^^ = 0.4 ~ 0.5fm, while the electromagnetic (EM) 

r.m.s. radius of the proton is \J < r"^ >p = O.Sfm. A formulation© to derive the 
short range part of the A^-A^ interaction directly from the quark-quark interaction 
has also been given. 

Now it seems to us quite natural that a possible physical origin of our repulsive 
core in the vrvr iso-singlet S'-wave state is the same mechanism as that in the nuclear 



force. Actually an interesting work from this viewpoint has been done by Ref.51). 
The values of our hard core radii in the tttt interaction are collected in Table IV with 
EM r.m.s. radius of the pion, in comparison with the mentioned cases of the a-a 
and the A^-A^ interactions. 



*•* This possibility had been noticed in Ref.E 
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Table V. Scattering length obtained from our fitting, in comparison with other estimations. 



Our estimation 

O-len—O.'len'^t: 


0.69fm (rc=0.68fm) 
0.32fm {r^=lfi^pn) 


Other estimation 


0.37±0.07fmp9 
0.34±0.13fm 


)|53|) 
54]) 


Current Algebra 
(Non-linear a model) 


3 

=0.22fmE3''E3) 



Hard core radii in the respective cases may be directly related with and may 
have similar orders of magnitude to the structural sizes of corresponding composite 
particles, considering our relevant mechanism. This expectation seems to be realized, 
as is seen in Table IV. 

However, one may feel that r'^^ seems to be rather large. We have made some 
analysis on this point. To reproduce the characteristic phase shift below KK thresh- 
old, a certain magnitude of the core radius is required. A permissible minimum value 
of the radius may be 2.5GeV~^. 

Here it is noted that our analysis wit h t he "soft-core" background Eq.(3-2) has 
been stimulated by the corresponding onel^ in the case of nuclear force. 

We have also estimated the vrvr scattering length aien=a,ff,n-'''c using the proper- 
ties of (J in Table I and the values of Tc from both in Tables I and HI. (We have 
used the form of relativistic B.W. formula, Eq.(2-14), in order to estimate a^^.) The 
results are given in Table V in comparison with the other estimations, which seem 
to be mutually consistent. However, we note that the values of a^^ depend largely 
upon the form of the relativistic B.W. formula, and effects due to various possible 
effects of hadron interactions (see §3.4.). 

(D) Possible existence of chiralons; new scalar and axial-vector meson nonets 

Classification of mesons usually follows the L^-coupling scheme of the non- 
relativistic quark model. We generalized it covariantly, keeping good L and S quan- 
tum numbers, to the boosted L5'-coupling scheme. 0^ In this scheme, the Bargmann- 
Wigner spinor wave function, which is a covariant generalization of Pauli spin func- 
tions, implies that the constituent quark and anti-quark are in "parton-like" parallel 
motion Ba), i.e., they move with the same velocity as that of the meson. "Usual" 
scalar and axial meson nonets can be classified in the first L-excited state in this 
scheme. 

On the other hand, in most ENJL models,© as a low energy effective theory of 
QCD, there exist four chiral symmetric nonets with J^'-''= 0"^"^, ^, 1 and 
In ENJL essentially only local composite quark and anti-quark operators are treated, 
thus missing L-excited states in principle. Here it is to be noted that the constituent 
quark and antiquark inside of these 0^^ and 1^^" nonets are shown mutually in 
an anti-parallel motion. Members of these nonets are apparently different from the 
above-mentioned Bargmann-Wigner particles 

Accordingly, it seems to us that, in order to complete the meson spectroscopy, 
we must take these extra "ultra-relativistic" 0"'"''' and l"*"*" nonets into account, in 



20 



S. Ishida et al. 



addition to the "non-mlativistic" ones described by the LS couphng scheme. We 
caU these "chiralons" . E^) 

It is expected that production of chiralons will be enhanced in pp central col- 
lisions, where constituent valence quarks and sea-antiquarks in anti-parallel motion 
exist. Thus if we assume that a is a member of scalar chiralons, enhancement of a 
in central collisions, mentioned in (B), may be clearly explained. 

The /o(1300) given in our analysis in §3 may be assigned to a member of the 
normal O"^"*" nonet in the ^Pq qq state. Also we are examining the possibility of 
/o(980) being a hybrid meson with a massive constituent gluon.® 

Finally we would like to note that one of the authors has examined 0) recently 
the properties of (t(555), obtained in this work, from the viewpoint of the effective 
chiral Lagrangian of the linear cj-model, and shown that C7(555) may be identified 
with the chiral partner to the Nambu-Goldstone vr-meson. 

§5. Summary 

In this paper we have re-analyzed the phase-shift between the vrvr-threshold and 
ISOOMeV by applying a newly developed lA-method. The results are summarized 
as follows: 

(1) The TTvr — > tttt phase shift has been analyzed in the mass region of the tttt- to 
the KK- thresholds, giving strong evidence for the existence of a, with mass of 
553.3 ± 0.5siMeV and width of 242.6 ± 1.2stMeY. 

(2) Our best fit value of a radius of the repulsive core in tttt interactions, whose 
introduction is the keystone for leading to the existence of (t(555), is 3.5GeV~^ (in 
case of "hard core" )/5.4GeV-^ ("soft core"). This value is almost equal to the 
structural size of the pion (3.5GeV~^). This situation is quite similar in the case of 
N-N interactions, giving support to our interpretation of the repulsive core. 

(3) The properties of /o(980) have been investigated from the data including those 
over the KK threshold. The mass is obtained as 993.2 it 6.5st i 6.9sysMeV. Its 
total width varies due to the uncertainty of KK coupling, to be 78.3 ± 11.5stMeV 
(fit A), 114.5 ± l5.5stMeV (fit B) or 344MeV (fit C), corresponding to the three 
different treatments of the elasticity and the vrvr — > KK phase shift, both of which 
have large experimental uncertainties. However, fit C has difficulty to reproduce the 
experimental vrvr KK phase behavior. 
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